Thermoelectric transport with electron-phonon coupling and electron-electron 

interaction in molecular junctions 
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Within the framework of nonequilibrium Green's functions, we investigate the thermoelectric 
transport in a single molecular junction with electron-phonon and electron-electron interactions. By 
transforming into a displaced phonon basis, we are able to deal with these interactions nonpertur- 
batively. Then, by invoking the weak tunneling limit, we are able to calculate the thermoelectricity. 
Results show that at low temperatures, resonances of the thermoelectric figure of merit, ZT, occur 
around the sides of resonances of electronic conductance but drop dramatically to zero at exactly 
these resonant points. We find ZT can be enhanced by increasing electron-phonon coupling and 
Coulomb repulsion, and an optimal enhancement is obtained when these two interactions are com- 
peting. Our results indicate a great potential for single molecular junctions as good thermoelectric 
devices over a wide range of temperatures. 

PACS numbers: 72.15.Jf, 72.10.Di, 73.23.Hk, 85.65.+h 



I. INTRODUCTION 

Recently, the potential afforded by nanoscale engineer- 
ing 1 has revitalized interest in developing novel thermo- 
electric materials for the generation and harvesting of 
energy. It is well accepted that nanoscale materials en- 
gineering in principle creates unlimited opportunities for 
the creation of more efficient energy-conversion devices 2 
and thus expands the potential of using thermoelectricity 
for meeting the challenge of being a sustainable energy 
source. 3 However, questions remain about the best ways 
for manipulating the microscopic properties of the mate- 
rial so that enhanced performance occurs. 3 ' 4 

The thermoelectric performance is typically character- 
ized by the figure of merit, ZT, 4,5 which is defined as 
ZT = G e S 2 T / k, where G e is the electronic conductance, 
S is the thcrmopower, T is the temperature, and k is 
the thermal conductance. Increasing the value of ZT in- 
creases the efficiency of heat-electricity conversion. The 
dependence of the figure of merit on both charge and 
energy transport shows that thermoelectric efficiency is 
strongly affected by the underlying electronic and vibra- 
tional properties of a material. These dependencies are 
especially transparent in molecular junctions, 6 as charge 
accumulation on the junction causes Coulomb interac- 
tions (e-e) to perturb the electronic structure 7 and the 
electron-phonon (e-ph) coupling to perturb the vibra- 
tional modes and the conformation of the junction. 8 

The energy scale of e-e interaction is usually much 
larger than that of e-ph interaction. However, at the 
atomic and molecular levels, the electrodes can screen 
the Coulomb repulsion, reducing it to the same order of 
magnitude as the e-ph interaction. The interactions now 
compete. Thus, it is of fundamental and practical im- 
portance to explore the effects of this competition to gain 



insights into the optimization of thermoelectric transport 
for better energy-conversion devices. 

In this work, we investigate the thermoelectric trans- 
port in a single molecular junction with e-e interaction 
and e-ph coupling of arbitrary strengths in the frame- 
work of the nonequilibrium Green's function (NEGF) 
method. 9 Although there are considerable efforts toward 
understanding the effects of the c-c 10 ~ 13 or e-ph 14,15 in- 
teractions, what happens to the thermoelectric transport 
when e-ph coupling competes with e-e interaction in a 
full range of strength is still an open question. In con- 
trast to previous work in the literature, 16,17 we treat the 
e-e and e-ph interactions within the molecular junction 
nonperturbatively, by a transform of the phonon basis 
with effective displacements. This treatment moves our 
use of the NEGF framework beyond the weak e-ph cou- 
pling perturbative analysis, 14,18 ' 19 the strong e-ph cou- 
pling limit of canonical transformation (the Lang-Firsov 
approach), 20-22 and the mean-field approximation in the 
strong Coulomb repulsion regime. 12,22 



II. METHOD AND APPROXIMATION 

We start with the standard Anderson-Holstein Hamil- 
tonian: 23,24 H = H mo \ + ifp + -Pleads- H mo \ describes the 
molecular junction of one orbital level, with Coulomb re- 
pulsion between electrons of opposite spin orientations 
and additional coupling to the vibration of itself, which 
is conventionally assumed to be: 9,16, 17,21-24 

a 
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In the first term on the right-hand side, and a create 
and annihilate a phonon with energy ujq while in the sec- 
ond term d£. and d a create and annihilate an electron of 
spin a at the molecular level with energies e-j- = Eq + Ae/2 
or £j_ = Eq — Ae/2. The third and fourth terms describe 
the e-ph and e-e interactions with strengths A CT and U. 
The expression 



H 



leads 



k,(7,u— L,R 



describes the left and right electrode leads with N£ a = 

^kau^kav the number operator for electrons of reservoir v 
with wave number k and spin a. The cxprcssoin 



a^kav 



k,(T,u— L.R 



describes the tunneling Hamiltonian of the electron hop- 
ping between the molecule and the electrode leads. Here, 
due to the large mismatch of vibrational spectra between 
the molecule and metallic electrodes, the phonon trans- 
port is not considered. 

Computing the transport for this model within the 
NEGF formalism requires the knowledge of various 
Green's functions for the different parts of H. We start 
by analytically solving the eigenproblem for the molecu- 
lar part of the junction. The Hilbert space of the molec- 
ular part is spanned by the basis n)}, where are 
the four possible electron states |0), \a), \a), \aa) and 
\n) denotes phonon states with n = 0, 1, . . . , oo. To non- 
perturbatively treat the e-ph and e-e interactions with 
arbitrary strengths, we first block diagonalizc the Hamil- 
tonian with respect to electron states: 



(0|# mo l|0) 
(cr|iJ mo i|cr) 

(aa\H mo i\<ra) 



wod^d, 



>a+ e a + A cr (d t + d), 



wo a 

ojoa^a + E £ <y + E + a) + U. 



|^moi|0) with conven- 
\n) = [(o+) n /Vni]|0). 



It is easy to further diagonalize 
tional Fock phonon states: |n) 
However, in order to diagonalize the other two matrix el- 
ements, we need introduce a new phonon basis, 25,26 with 
displacements shifted by different electron states through 
the e-ph coupling: 



l">« = UlT/Vri] cxp (- 5 2 /2 - g*tf)\0), (1) 



where 



it. 



- denotes the new creator that creates 
a phonon displaced from the original position by a value 
g$ depending on the electronic state, that is, g = 0, 
g <y = Acr/wo, and g aS = g a + g s . Clearly, when electrons 
are absent on the molecular quantum dot, g = and the 
displaced phonon basis reduces to the normal Fock state 
of the phonons. We then can write 

H-ffmolIu) = UJoA^Aa- + e a - wogl, 



(aa\H mo i\a<7 



(7(7 



+ e a + e s - woflW + U. 



Therefore, with the help of the new phonon basis, the 
solution to the eigenvalue problem is 

-£fmoi|0,"-) = nLL> Q \0,n} , (2) 
H mol \<7,n) a = (nu>Q + i a )\(j,n) a , (3) 
H mo \\<ra,n) ai j = {ntoa + e aS + U)\aa, n) aS , (4) 



where 



e a - uj g o 



2u g a g s . 



The negative term —2ojQg a g s evidences the attractive in- 
teraction between different electron states induced by the 
e-ph coupling. 

With -ff mo i diagonalized, we can now analytically cal- 
culate the advanced and retarded Green's functions of 
the molecule, which are found to be (see Appendix A for 
details) 



z 



n,m— 



o — f3mu> 



(nui +i aa +U) 



DIM, (5) 



where 



A ™ = (n-m)wo +£a, 

= (n - mjojQ + (e<j - 2ui Q g a g s + U), 
D nm {g a ) = (-l) m (m|n) CT = {-l) m s (m\n) a s 

min{n,m} 



= e -9l/2 J2 
fe=0 



(-l) k yM^-ga +m ~ 2k 
(n - fc)!(m - k)lk\ ' 



Z = (l + N ph )(l + e~ ps ' + e~ pe ' + e -^ +u ^). 

Here N p h = l/(e^ LJ " — 1) denotes the Bose distribu- 
tion of the phonon population with inverse temperature 
/3 = ksT. Note, in the zero-temperature limit, Eq. (5) 
is consistent with the "atomic limit" Green's functions 
in Rcfs. 24,27 through using canonical transformation. 
The advantage of our method is that at finite temper- 
atures our results still give Green's functions analyti- 
cally and explicitly, while for the atomic limit with the 
Lang-Firsov canonical transformation, the expressions of 
Green's functions contain the electronic level occupa- 
tion, the value of which must be found through a self- 
consistency iteration. 

With the molecular part treated nonperturbativcly, we 
now follow standard paths to build the Green's func- 
tions of the molecule-lead-electrode system. By using the 
Dyson equation and the Keldysh formula, 9,21 we have the 
total retarded (advanced) Green's function 

r< r i a ) — \tr<r{a)\-i v r ( a ) l- 1 
and the total lesser (greater) Green's function 



G 



<(>) 



-.<(», 

'tot,cr i -'tot,o- ' 
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The total self-energy has two contributions: Stot.c = 
^lead a + ^int a i where S; nt a is the contribution from 



the e-ph and e-e interactions, following from 9,21 £ 



r(a) 



G 



r(a)~ 1 _ r ,r(a)- 1 
0,cr ) 



andE^J = G r a - l G^ {>) G a a -\ where 



Gq " (w) = (oj — e a ± i0 + ) _1 denotes the noninteracting 
Green's functions without involving e-ph and e-e inter- 
actions. Here £i ca d,er depicts the contribution from two 
tunneling parts between the molecular quantum dot and 
leads. 

So far, no approximations have been made for the 
representation. However, to exactly obtain Ei oa d,cr hi a 
strong correlated system is highly nontrivial. In the fol- 
lowing, we take the lowest order approximation of 5]iead, a 
as in a noninteracting system and obtain 9 



lcad,(T 

^lead.tr 

y,r(a) 
lcad,<r 



-i\rZ(i-f L ) + r?(i -/«)], 



where T a = (T L a + T*)/2, T%(u>) = 2tt J2 k |CI 2 <K" - 
e ka .), denoting the molecule-electrode coupling functions, 
are energy independent in the wideband limit, and /„ = 
.\-\y- 1 are th e Fermi-Dirac distributions of two 
reservoirs. This approximation is valid in the weak tun- 
neling limit; i.e., we take small values of Y v al comparing to 
all other energy scales. This weak tunneling approxima- 
tion is consistent with the polaron tunneling approxima- 
tion in Ref. 28, under the Lang-Firsov canonical trans- 
formation. Note this weak tunneling approximation also 
ignores the Kondo effect, which is justified since we are 
interested in regimes beyond the Kondo temperature. An 
interpolative approach for self-energies 27 may be used to 
relax this weak tunneling limit. 

With the Green's functions in hand, we can now study 
the quantum transport by calculating the charge current 

J v = -e(J^N^) 

k<7 

and heat current 

leaving electrode v, which in terms of the total Green's 
functions are 9 

Jv = ~jT E J ^ C^^tot.CT + fv^[G[ot,o ~ ^tot.o-]) ! 

r Q = I E / dtu ^ ^ K G fot ia + mez*,, - Gt t, 

After some algebra (see details in Appendix B), we find 
that the currents going from the left electrode to the 
central system are 



where 



with 



7» =7}M+7L(w) 



L r Ji. 



r^r 



(8) 



We can also obtain the same expressions for Jr and Iq 
with L <-> R, depicting the electron and heat current 
from the junction to the right electrode. 

Any meaningful transport theory in terms of NEGF 
must respect current conservation: 9 



e f dtu m r 



int^tot 



0. 



The self-consistent Born approximation is an exam- 
ple of a conserving approximation. Without the self- 
consistency of Green's functions and self-energies, the 
current nonconserving issue Jl =/= —Jr generally ex- 
ists, for example, in the perturbation approximation of 
weak e-ph coupling 29 and the canonical transformation 
for strong e-ph coupling. 20 People usually choose to cal- 
culate the symmetrized current J = (Jl — Jr)/2 to avoid 
the current nonconserving issue. The present approach, 
however, satisfies Jl = — Jr directly. It also avoids the 
time-consuming self-consistent refinement in the mean- 
field treatment of Coulomb repulsion. 9 

Energy conservation also gives the relation — cJo-^ph = 
Iq + Iq- 14 For linear response, the system is in equilib- 
rium such that A^ph = and Iq = —Iq, since N p ^ follows 
the Bose distribution for given temperature T. There- 
fore, the equilibrium temperature is directly used for cal- 
culation, and indeed our heat currents satisfy Iq = —Iq- 



III. THERMOELECTRIC TRANSPORT AND 
DISCUSSIONS 

The thermoelectric coefficients are conventionally con- 
sidered around the linear response region: 21 



Hl(r) = A*f ± eAV/2, 
T L{R) =T± AT/2, 



which yields 



J) 



J \ _ ( G e G e ST \(AV\ 

Iq)~\ G e ST (G e S 2 T + k e )T J { AT/T J 



where 



G e = J/AV\ AT=0 = e 2 C 



(9) 



Jl = \J ^T(u>)[f L (u,)-f R (u>)], 



(6) 



with 



1 f du> 



Iq = I I - Ml)TH[/ l H - f R (w)], (7) 



S = -AV/AT\ J=0 = £i/(er£o), (10) 
Ke = Iq/AT\ j=0 = (C 2 - C 2 /£ )/T, (11) 

C n = ~ J dwT(w)(w - vL F ) n . (12) 
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FIG. 1: (Color online) Electronic and thermal conductances 
and Lorenz ratio as a function of fiF for different e-ph and 
e-e interaction strengths: so = —2hwo, Ae = 0, = = 
0.05hujo, huo — 30 meV, and T = 35 K, which are typical 
experimental values. 18 



While G e S 2 T + K e is the thermal conductance at zero 
voltage bias, n e is the conventional thermal conductance 
of electrons at zero electron current. Another important 
quantity is the Lorenz ratio L = K e /(G e T). For macro- 
scopic conductors, the Wiedemann-Franz (WF) law re- 
lates the electronic and heat conductances via the uni- 
versal relation L = L = (7r 2 /3)(fcs/e) 2 , which indicates 
that charge and energy currents suffer from the same 
scattering mechanisms such that more electrons carry 
more heat and vice versa. In a single-electron transistor, 
the Coulomb blockade effect leads to the strong violation 
of WF law: L/Lq ^s> l. 30 However, in order to obtain a 
large ZT = G e S 2 T/n e = S 2 /L, the opposite violation of 
WF law, L/Lq <C 1, is desirable. 

Figures 1(a) and 1(b) show the electronic conductance 
G e as a function of fXp for different values of A and U. 
We see that the number of resonance peaks increases 
as the e-ph coupling increases: Large A excites more 
phonons and enables multi-phonon-assisted tunneling. 
At low temperature, we can approximate G e as T{^p) 
since —df/duj as 8{uj — [ip). Therefore, from Eqs. (5) 
and (3), the positions of the resonance peaks are de- 
termined by the poles of Green's functions; i.e., /ip of 
resonance positions are equal to A^m or A n 2 m . How- 
ever, exponential functions and D nm in Eq. (5) weight 
the resonances and make the peaks unobservable for cer- 
tain parameter ranges. From Eq. (5), when U is so weak 

(2) 

that e aS + U < e a < 0, A„„ < thus dominates the 
resonance peaks [see Fig. 1(a)]. While U increases up to 
£<?o + U > s (j , two resonance branches, A nm > and 
Anm < 0, emerge [see Fig. 1(b)]. For large t/, when 
A increases further such that again e aS + U < e a , then 



A nm < redominates the resonances and two resonant 
branches merge together, as shown in Fig. 1(b). It is 
a consequence of the competition between the e-ph cou- 
pling and e-e interaction with the molecular quantum dot 
junction. 

In Figs. 1(a) and 1(b), we also observe that increasing 
e-ph coupling decreases the peak value of the main reso- 
nance of G e but increases the values of side peaks of G e . 
This occurs because D 2 = e~ 9 decreases with g increas- 
ing, while for a positive integer n, D 2 n = g 2n e~ 9 /n! has 
the opposite tendency at g G [0,-y/n]. When e-ph cou- 
pling g increases further (g > \/n), D^ n decreases again 
so that side peak values of G e will be repressed by the 
strong e-ph scattering. Comparing Figs. 1(a) and 1(b), 
we further see that for the weak and moderate A, increas- 
ing the Coulomb repulsion reduces G e , which is a conse- 
quence of the factor e~@ u in Eq. (5). While for strong 
e-ph coupling, the Coulomb repulsion mainly shifts the 
positions of the spectrum of G e while leaving its magni- 
tude almost unaffected. This is because for large A, while 

(2) 

the resonance positions A nm depend on U, the Green's 
functions are dominated by D nm , which, however, is U 
independent. 

The A and U dependence of thermal conductance n e is 
similar to that of G e , as shown in Figs. 1(c) and 1(d). We 
note that the resonance positions of n e do not coincide 
with those of G e but instead coincide with the valleys 
of G e . This arrangement reflects the different ways in 
which the inelastic scattering induced by e-ph coupling 
and e-e interaction degrade heat and electrical currents. 
In fact, around the resonances of G e , from Eq. (12) it is 
clear that — 0, so we have small values of n e from 
the definition Eq. (11). As a consequence, we obtain 
the strong violation of WF law L/Lq <C 1 around the 
resonance points of G e [see Figs. 1(e) and 1(f)]. 

Because ZT = S 2 /L becomes large as the Lorenz ra- 
tio goes to zero while the thermopower remains finite, we 
expect large ZT peaks will emerge around the resonance 
points of G e . However, as shown in Fig. 2, large ZT oc- 
curs only at the sides of the resonances and drops back 
to zero dramatically at exact resonance positions. This 
occurs because the particle-hole symmetry at the reso- 
nances [Ci ~ from the definition Eq. (12)] zeros the 
thermopower [S ~ from the definition Eq. (10)]. More- 
over, increasing A increases ZT as well as the number of 
ZT peaks, which shows that large ZT is favored by multi- 
phonon-assisted tunneling. In addition, the Coulomb re- 
pulsion increases ZT such that optimal ZT is obtained 
at the merging regime of two resonant branches, which 
is a consequence of the competition of e-ph coupling and 
Coulomb repulsion. In other words, the optimal ZT is lo- 
cated at £<rcr + U = e a as we discussed above in Figs. 1(a) 
and 1(b). Taking e a — e s = —2hu>Q,U = AHljq as in 
Fig. 1(b), we then can predict that the optimal ZT can 
occur when we choose A CT = y/2/3hu>Q tn O.SlhujQ, which 
is indeed the case shown in Fig. 2(d). In turn, given the 
molecular level energy and the e-ph coupling strength, 
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FIG. 2: (Color online) S and ZT as a function of fj,p for 
different A and U. The Hamiltonian parameters are the same 
as those in Fig. 1. 
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FIG. 3: (Color online) Temperature-dependence of ZT for 
moderate A = 0.5&jo, with (a) U = and (b) U = 4fiwo- 
Other parameters are the same as those in Fig. 1. 



tune the conduction electron density. Additionally, there 
is also the flexibility of adjusting the coupling geometry 
between molecule and leads. 



IV. CONCLUSIONS 

In conclusion, by a transformation of the phonon basis 
we were able to nonperturbatively deal with the molec- 
ular quantum dot for arbitrary e-ph coupling and e-e 
interaction strengths. After analytically calculating its 
Green's functions, we coupled the molecular quantum 
dot to the electrode leads in the weak tunneling limit, 
and then computed the thermoelectric transport prop- 
erties numerically. We studied the synergistic effect of 
e-ph and e-e interactions and showed that at low tem- 
peratures large ZT occurs at the sides of resonances in 
electronic conductance but drops dramatically to zero at 
resonance. We found that increasing e-ph and e-e in- 
teractions increases ZT, although with G e repressed. In 
particular, large ZT is favored by multi-phonon-assistcd 
tunneling. More interestingly, we found that an optimal 
ZT emerges when these two interactions were competing. 
Finally, we showed that a large ZT can be obtained in a 
wide range of temperatures. 

It would be interesting to consider the Zeeman split- 
ting of the molecular orbital energies produced, for ex- 
ample, by an external magnetic field or ferromagnetic 
leads. Spin-related thermoelectric effects may aid the 
optimal design of novel thermal-spintronic devices and 
singlc-molecule-magnet junctions. 31 Extending our non- 
perturbative approach to electronic coupling to multiple 
vibrational modes 33 will be an interesting topic. It would 
also be desirable to combine the present method with 
some ab initio electron structure theory, like the den- 
sity functional theory within the local density approxima- 
tions, for more realistic calculations. 9 Finally, we would 
like to remark that calculating the exact self-energy Si ea d 
remains an important open question, which deserves in- 
vestigation in the future. 



we also can choose a proper e-e repulsion strength U to 
optimize the efficiency of thermoelectricity. 

Finally, the temperature dependence of ZT is illus- 
trated in Fig. 3. Interestingly, we see that increasing T 
first decreases and then increases ZT. Coulomb repul- 
sion enhances ZT at low temperature but suppresses it 
at room temperature. Nevertheless, ZT remains large. 
This result indicates a great potential for single molecu- 
lar junctions as good thermoelectric devices over a wide 
range of temperatures. 

In view of the substantial progress that has been made 
in the area of molecular devices, 32 the findings in the 
present work will open a new avenue for using molecular 
junctions for the optimal design of efficient energy conver- 
sion devices. To experimentally observe the predicted ZT 
enhancement, the molecular junction should be gated to 
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Appendix A: Calculation details of nonequilibrium 
Green's functions 

Before calculating the retarded (advanced) Green's 
function, we first detail the calculation of the lesser 
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Green's function G<(t) = i(dl(0)d a {t)) in the frequency 
domain: 

r+co 

G<(cu)=i / dte^{dl(Q)d a (i)) 



+ 00 



-PE V -i(E v -E^)t 



x (^4(0)1^)^1^(0)1^) 



PE V 



If 1p 



x ^|4(om(^k(oM, 



(Al) 



where we used p = e~P H -™ l jZ with Z = Tr(e"' 3ff '"° 1 ) 
: i + N pb )(l + e-' 3i " + e"^" 5 + e -^°*+u)) an d 7V ph is 
the Bose distribution N p h = l/(e l3Ld0 — 1). Here |<p) and 
\ip) are the possible eigenstates |0,n) , |er, n)^, |ct, n)g-, 
I era-, n) aS , and i?^ and £/,/, are the corresponding possible 
eigenvalues uluq, tiluq + e a , uujq + uluq + E aS + U. 
There are only two nonzero combinations of \(p) and 
for calculating G<(u>): (1) \<p) = \<J,n) a and = 
|0,m) , or (2) \tp) = \aa, n) a9 and = \a,m}&, such 
that the lesser Green's function can be reduced to 



00 00 



x e-^" uo+ ^^(n|m) (m|n) cr 

+ <5(w - (n - m)w - (e^ - £ff + £/)) 

x e-^( n " 0+ ^ s+c/ ) aff ( ? i|m) ff *<m|n) ffff ]. (A2) 

The detailed expression of b(n\m) Cl denoting the inner 
product of modified phonon states with effective displace- 
ments gb and g c , can be derived as follows: 

, (a + g b ) n 



b (n\m) c = (0|- 



■ exp (-g b /2 - g b d) 



exp(-g c 2 /2- : g c a))\Q) 



exp [-(ff b - ge ) 2 /2] 
\/n!rn! 

x (0|(a + ffb )™e(- 9cat )e ( " 9i ' a) (a t + <? c ) m |0) 
cxp[-(.9fa - ,g c ) 2 /2] 



V n\m\ 

x(0\(a + g b -g c r(tf+g c -g b ) m \0) 
exp[-(g b -g c ) 2 /2] 

min{n,m} 

x £ k\C*{g b -g c ) n - k C* n ( 3c - !n ; i 

fc=0 

(-ir^ nm (. 9b -.g c ) 



rn — k 



(A3) 



where 



i{n,m} 



771 — 2 fc 



fe=0 



(-l) fc \/n!TO!a;™ + 
(n-/c)!(m-fc)!fc! 



is invariant under the exchange of indices n,m. Note, 
to get the third equivalence, we utilized the relation 
exp {cd)f(d\ a) = f(a) + c, a) exp (ca). 

Therefore, the lesser Green's function can be further 
reduced: 



2m 



OO 



A( 2 2je-^+ E - +c/ >l £ 2 m (g CT ), (A4) 



where 



A im = {n-m)u +e*, (A5) 
ASS, - {n-m)uj + (i a -2uj g a g s + U). (A6) 

Similarly, for the greater Green's function G>(t) = 
—i(d a (t)dl.(0)), we can obtain 



G>(u) = -i J dte wt (d„(t)4(0)) 

= -i f dte'-*^^(^|p e ^"- t d ff (0)e-^' t |^ 

x(v>|4(o)M 

x(<p\d<rW}(lP\dU<p) 

= -^E[^- A ^) e ~" mwo + 

5( W - A« ) e -^"-o+e- ff )] £) 2 m (g CT ). (A7) 

Note here the two nonzero combinations of \ip) and for 
calculating G> (a; ) are (1) \ip) — |0,m) and — |a, n)^, 
or (2) \tp) = \a,m)a and = \aa,n) aS . 

Then, following the relation G r {t) = Q(t)(G> (t) - 
G<(t)), G a {t) = -e(-t){G>(t) - G<(t)), and utilizing 



0(i) 



27ri cj — ?'0 + ' 

we have the retarded (advanced) Green's function: 



±1 f^l [ d t e ^^—^i 
2tt J 2iri J lu 2 T «0+ 



x[G>(wi)-G<(wi)] 
27ri wi — =p i0 + 



(A8) 



Substituting the expressions of the greater and lesser 
Green's functions, we arrive at Eq. (5). 
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Appendix B: Derivations for the current expression Substituting the expressions of various noncquilibrium 

Green's functions, 

Here we detail the calculation of the current through 

the interacting system. The electronic current from ^ 

left contact to central system is defined as J L = ^tota = — ^ > (B2) 

— e (Sfctr dN^ a /dt), which is generally reexpressed as 9 {G a )~ 1 ± i^a 

/"<<(>) f~lr y<(>) f a tTZ1\ 

„■„ r A, , ^tot.o- ~~ "tot.o-^tot.cr^tot.o-) \ DO ) 



Jl — J — / y [r^tot.cr + r^/i(G^ ot)CT - G" ot i0 . ;j • 

(Bl) we have 



< 

tot, a 



2tt ^ 



*w) - (2*r CT + (G^)- 1 - (Gj)- 1 )^^ 



do; x - 



r^< 



2 7r ^[(G;)- 1 +zr (T ][(G-)- 



ie f d^^T^G< 



r^/ £ (G^-Gj) 



^ 27r^[(G;)-i + ir (r ][(G«)-i-ir <r ] 



2tt ^ [1 

a L 



zr^G^li-zr^Gj] 



J 



V L J L {G> a - G<) 



duj v r^Gg 

27r z^ [l + ir^li-ir^G-] 



(B4) 
(B5) 
(B6) 



(B7) 



In the integration of the last line, the lesser and greater 
Green's functions contain the Dirac i5 functions, which 
only have finite nonzero values at the resonant points. 
However, at those resonant points, the retarded and ad- 
vanced Green's functions have divergent values, which 
finally lead to zero integration values. Therefore, the 
contribution of the last integration is zero, and we finally 
arrive at the electron current expression, Eq. (6). It is 
easy to get the same expression for Jr with R <-> L, such 
that current conservation is explicitly preserved. 



Following the first law of thermodynamics dQ = dE — 
/j,dN, we have the current relation: Iq = Q = E — fiJ. 
Therefore, following the similar calculation, the heat cur- 
rent is straightforwardly obtained as Eq. (7). We can ob- 
tain the similar expression for the heat current through 
the right reservoir Iq. Based on the electron current ex- 
pression and the expression of heat current carried by the 
electron, we are capable of investigating the thermoelec- 
tric transport properties. 



* Present address: Theoretical Division, Los Alamos Na- 
tional Laboratory, Los Alamos, New Mexico 87545, USA. 
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